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ABSTRACT 

We present first results of a spectroscopic survey targeting iGselected galaxies at z = 2.0 — 2.7 
using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength 
coverage of 1.0-2.5 /mi for 26 if -bright galaxies (K < 19.7) selected from the MUSYC survey using 
photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest- 
frame optical emission lines or the redshifted Balmer/4000A break. Twenty galaxies have spectroscopic 
redshifts in the range 2.0 < z < 2.7, for which bright emission lines like Ha and [O in] fall in 
atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The 
median 2a upper limit on the rest-frame equivalent width of Ha for these nine galaxies is ^10 A. 
The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population 
models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with 
zero within la for the remaining four. Thus, both the Ha measurements and the independent stellar 
continuum modeling imply that 45% of our iGselected galaxies are not forming stars intensely. This 
high fraction of galaxies without detected line emission and low SFRs may imply that the suppression 
of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy 
formation models. However, obscured star formation may have been missed, and deep mid-infrared 
imaging is needed to clarify this situation. 

Subject headings: galaxies: high redshift — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Observations imply that stellar populations of high- 
mass galaxies were formed at hi gher redshift than 
those of low-mass galaxies (e.g., ICowie et all 119961 
Uuneau et al J 12005). Recent hierarchical CDM models 
are able to produce dead, massive galaxies at high red- 
shift by incorpor ating feed-back from active galactic nu- 
clei I' AON: e.g iCroton et all l2006f iBower et all 120061: 
iDe Lucia et al.ll2006t iHopkins et al.ll2006j) . In order to 
observationally determine when and how star formation 
in massive galaxies was suppressed it is necessary to iden- 
tify and study these objects out to the highest redshifts. 

Recently, massive and apparently dead galaxies have 
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been identified at z > 1.5 (e.g iMcCarthv et alJ | 2004 
[ Saracco et alll2005tlD~addi et al.l2005tlLabbe et aljhoosl 
iReddv et al.ll2005l 120061 iPapovich et al.ll2006jL Bevond 
z = 2 most studies rely on photometric redshifts and 
broad-band colors to identify these galaxies. However, 
as dust and age have similar effects on the broad-band 
spectral energy distribution (SED), the star formation 
rates (SFRs) are often not well constrained. Further- 
more, their photometric redshifts are not well calibrated, 
as only a few at z > 2 have spectroscopic redshifts 
ijDaddi et al J 12005). Thus spectroscopic redshifts and 
independent stellar population diagnostics are needed to 
determine the prevalence of "red and dead" galaxies be- 
yond z = 2. 

Spectroscopic confirmation of non-star-forming galax- 
ies at z > 2 is complicated due to their faint rest-frame 
UV emission and lack of nebular emission lines. Deep 
near-infrared (NIR) spectroscopy provides the best op- 
tion to confirm such galaxies at z > 2, as their rela- 
tively bright rest-frame optical luminosity allows direct 
detection of the stellar continuum. The optical contin- 
uum shape, and in particular the Balmer/4000 A break, 
can be used to derive redshifts for galaxies without emis- 
sion lines, and p rovides independe nt constraints on stel- 
lar populations ijKriek et al.ll2006j) . 

To study a high-redshift spectroscopic sample that is 
not biased towards galaxies with bright emission lines, 
we are conducting a NIR spectroscopic survey of K- 
selected galaxies with photometric redshifts z ~ 2.3. 
Here we report on a surprising result of our survey: the 
large fraction of galaxies with no detected emission lines. 
Throughout the letter we assume a ACDM cosmology 
with n m = 0.3, n A = 0.7, and H = 70 km s" 1 Mpc" 1 . 
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Fig. 1. — NIR spectra (purple) and optical-to-NIR photometry (blue) of the nine galaxies at 2.0 < z < 2.7 for which we detected no Ha 
emission. The upper panels show the ID original spectrum. The 'low resolution' binned spectra (160 A per bin) and the photometry are 
presented in the lower panels. All fluxes are given in 10~ 19 ergs s _1 cm - 2 A -1 . Regions with low or variable atmospheric transmission 
or with strong sky line emission are indicated in light gray. The best fit to the optical photometry and low resolution spectrum, allowing 
all values for Ay, is overplotted in dark gray. The best-fit model parameters are printed in each panel. The corresponding errors are the 
68% confidence intervals, derived from 200 Monte Carlo simulations. The \ 2 values are given per degree of freedom (JVdeg = 28). Models 
without dust generally provide good fits as well (light gray), and imply ages that are about a factor of 2 higher than the listed ages. All 
galaxies without detected H« emission are best fitted by evolved stellar population models with low specific SFRs. 



All broadband magnitudes are given in the Vega-based 
photometric system. 

2. GALAXY SAMPLE AND DATA 

The galaxies presented in this letter are drawn from 
the Multi- wavelength Survey by Yale-Chile (MUSYC). 
This survey provides optical and deep NI R photometry 
for several southern and equatorial fields ijGawiser et al.l 
20061 iQuadri et alJ l2006). The targets were selected in 
K (K < 19.7) to reduce the dispersion in stellar mass 



and to ensure adequate S/N in the NIR spectra. Addi- 
tionally, we required a photometric redshift in the range 
2.0 < z < 2.7, for which bright rest-frame optical emis- 
sion lines such as [O in] and Ha fall in the H and K 
atmospheric windows. The photometric reds hifts are de- 
rived followi ng the procedure described in lRudnick et all 
(200ll l200l) . 

We observed 26 galaxies with the Gemini near-infrared 
spectrograph (GNIRS) in 2004 September (GS-2004B- 
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Fig. 2. — The equivalent width of Ha vs. the specific SFR derived 
from our model fits to the spectra. Filled red circles are galax- 
ies with no detected Ha emission and purple squares are emission 
line galaxie s in our sample. Blue crosses are UV-selected galax- 
ies by lErb et alj 12006albld) . DRGs are indicated by open green 
diamonds. Upper limits for Ha and the specific SFR are 2a. Ex- 
pec ted relations between Wnr* and specific SFR are derived from 
the IBruzual fc Charlotl 120031) and IKennicutj 119981) models and 
drawn for a TiooMyr (gray dotted line) and a TxoooMyr {gray dashed 
line) model for the first 3 Gyr. Both Ha measurements and the in- 
dependent SED modeling demonstrate that in 9 out of 20 galaxies 
in our sample the star formation has been strongly suppressed. 



Q-38), 2005 May (GS-2005A-Q-20), 2006 January (GS- 
2005B-C-12) and 2006 February (GS-2006A-C-6). The 
spectra of two of these galaxie s have alrea dy been pre- 
sented by Ivan Dokkum et aD l|2005|) and fariek et"afl 
(2006). All galaxies were observed in cross-dispersed 
mode, in combination with the short wavelength cam- 
era, the 32 line mm" 1 grating (R=1000) and the Of! 675 
by 6'.'2 slit. In this configuration we obtained a wave- 
length coverage of 1.0 - 2.5 /im. The galaxies were ob- 
served for 1-4 hours, depending on the brightness of the 
target and the weather conditions. The observational 
techniques and red uction of the GNIR S spectra are de- 
scribed in detail bv lKriek et alJ l)2006j) . For each galaxy 
we extract a one-dimensional original and low resolution 
binned spectrum. 

To derive the stellar population properties and obtain 
redshifts for galaxies without emission lines, we fit stel- 
lar population models to the low resolution continuum 
spectra together with t he UBVRIz fluxes , following the 
technique described in iKriek et alJ l)2006f) . We use the 
IBruzual fc Gharlotl ((2003) models with a set of exponen- 
tially declining star formation histories, a Sal peterl l)1955f) 
initial mass function (IMF) betwe en 0.1 and 100 Mp, , 
and solar metallicity, and adopt the lGalzetti et all ( 2000) 
reddening law. The assumed model parameters (IMF, 
reddenin g law, metallicity) are i dentica l to those used 
by e.g ., lForste£_Sch£gIber etajJ lj2004j) . IShaplev et al.l 
( 2005f, and lPapovich et all l|2006T) . 

We obtained spectroscopic redshifts for all 26 galax- 
ies using rest-frame optical emission lines or the 



Balmer/4000A break. The 'break' redshifts have a me- 
dian uncertainty of |Az|/(l + z) = 0.017, as determined 
from fitting the low-resolution continua of emission line 
galaxies with z as a free parameter. 20 of 26 galaxies 
have spectroscopic redshifts in the range 2.0 < z < 2.7, 
for which Ha falls in the if-band. In what follows we 
restrict the sample to the galaxies in this redshift range; 
the full sample will be described elsewhere (M. Kriek et 
al., in preparation). We note that the six galaxies that 
fall out of this redshift range have z = 1.75 — 1.95. 12 

3. SUPPRESSED STAR FORMATION 

Surprisingly, nine out of the 20 galaxies in the sample 
show no emission lines in their rest-frame optical spectra. 
The spectra and best fits of these galaxies are presented 
in Fig. ^ For these galaxies we derived upper limits 
on the Haequivalent width (Whq) as follows. We have 
drawn 200 random redshifts from the redshift probabil- 
ity distribution, and determined the 2a upper limit of 
Wn a in each case from the measured noise properties, 
assuming a rest-frame H a FWHM of 500 km s _1 (see 
Ivan Dokkum et aT]l2004D and the best-fit stellar contin- 
uum. The adopted limit is the maximum value found 
in the simulations, excluding the highest 5%. The me- 
dian Wn a upper limit of these nine galaxies (corrected 
for Balmer absorption) is 10 A. 

Wn a is a measure of the ratio of current to past star 
formation, and the limits on Wn a may imply very low 
SFRs in these galaxies. We investigate this in Fig. |21 in 
which we plot the rest-frame Wn a upper limits (corrected 
for Balmer absorption) versus the specific SFRs (SFR / 
M*) derived from our model fits to the spectra. Remark- 
ably, all nine galaxies without detected emission lines 
(filled red circles) are best fitted by stellar population 
models with low specific SFRs, and the data points are 
broadly consistent with the expected relations between 
these properties (jKennicuttl 119981 IBruzual fc GharloU 
2003). For five galaxies the best-fit SFR is zero and 
the four remaining galaxies have best-fit values that are 
consistent with zero within la. The midmean 13 specific 
and absolute SFRs are 0.004 Gyr" 1 and 0.9M s /yr re- 
spectively, significantly lower than the 0.56 Gyr" 1 and 
128Af Q /yr found for the eleven emission line galaxies. 
Thus, both the Ha measurements and the stellar contin- 
uum modeling imply that the star formation in these nine 
galaxies has been strongly suppressed. We note that we 
find a similar relation when we plot Lh q vs. the modeled 
absolute SFR, as our galaxies span only a small range in 
stellar mass (0.9 — 4.6 x 10 11 Mq). These nine galaxies 
have a median stellar mass of 2.6 x 10 11 M Q , a median 
J — K color of 2.45, and have undergone a median of 
21 age/r e-folding times. Six ou t of nine are distan t red 
galaxies (DRGs, J - K > 2.3, iFranx et al1l2003|) . As 
can be seen in Fig. |21 several of the emission line galaxies 
are also best fitted by stellar population models with low 
specific SFRs. This may suggest that the gas in these 
galaxies is not ionized by hot stars; we will explore this 
in a future paper. 

Formally, we find high best-fit values for the dust con- 
tent for eight out of nine galaxies without emission lines. 

12 The median of the differences between spectroscopic and pho- 
tometric redshifts (z spe c — 2 p hot)/ (1 + ^spec) is only —0.001 for the 
full sample of 26 galaxies. 

13 mean of the central two quadrants 
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However, Ay is poorly constrained for most of these 
galaxies, and models with zero or only small amounts of 
dust are consistent within la. We refitted all 9 galax- 
ies allowing only models without dust. The best fits 
(Fig. ^ light gray) also yield low specific SFRs, rang- 
ing from 0.001 — 0.013 Gyr -1 , with a median value of 
0.002 Gyr -1 . The median best-fit stellar age is 0.9 Gyr, 
which is a factor of ~ 2 larger than when allowing dust. 
While we cannot draw firm conclusions about the dust 
content in these galaxies, we note that high Ay may in- 
dicate that they are still in the process of losing their gas 
and dust. 

Figure El a l so shows the UV-selected z ~ 2 galaxies 
of lErb et alJ l)2006albl lcft. These galaxies have similar 
Wn a and specific SFRs as the if-selected galaxies with 
emission lines. However, there is no overlap with the 
if-selected galaxies without detected Ha emission; both 
the Who and the modeled specific SFR are higher for the 
UV-selected galaxies. Alt hough UV-selection is able to 
find m assive galaxies (e.g.. lShaplev et al.l2005l:lErb et al.l 
2006b), it does not sample the full distribut ion of their 
properties (see also Ivan Dokkum et al]l200^ . 

4. DISCUSSION 

We find that nine out of 20 if-selected galaxies at 2.0 < 
z spcc < 2.7 have no detected Ha emission (Wn a ^5 10 A) 
and are best fit by stellar population models with low 
specific SFRs (~ 0.004 Gyr -1 ), implying a fraction 
of galaxies with strongly suppressed star formation of 
45^2 The quoted uncertainty is derived assuming 
Poisson statistics, and does not include the following 
systematic errors and caveats. First, our z p hot selec- 
tion criterion could introduce biases, as systematic errors 
in photometric redshift may correlate with SED type 14 . 
Also, the if -band selection criterion could bias our sam- 
ple, as for starburst galaxies strong emission lines can 
contrib ute significantly to the if-band flux (Er b et ahl 
2006a). We do not expect this bias to be strong, as the 
median contribution of the emission lines to the if -band 
is only 0.04 mag for the eleven emission line galaxies. 
Third, incompleteness may play a role, as we observed 
only ^20% of the galaxies that meet the selection crite- 
ria. We note, however, that according to RS- and KS- 
tests our if -selected sample has a similar distribution of 
rest-frame U — U-colors as the large mass- limited sample 
(> 10 n M Q ) bv lvan Dokkum et al.l <|2006|) . when apply- 
ing the same if -magnitude cut. Furthermore, the as- 
sumption that Ha emission is only due to star formation 
may lead us to underestimate the fraction. 

Finally, and perhaps most importantly, we may miss 
star formation with very high rest-frame optical extinc- 
tion. Although most local ULIRGs have high integ rated 
Wn a (W Ha +[NH] ~ 90A. iLiu fe Kennicuttl Il995f) . for 
some the star-burst regions are almost completely ob- 
scured (e.g., Arp 220 has WHa+[Nii] = 18 A) and these ob- 
jects might be misinterpreted in our analysis. Available 
24 /jm imaging on DR G samples shows that 30-50% have 
no MIPS counterpart iTWebb et alJ l200fit IP apovich et al.l 
120011 iReddv et al.ll2006ft . and Redd v et alJ (|2006ft finds 
that red galaxies with low MIPS fluxes typically have 

14 This effect may be small, as the distribution of specific SFRs 
of the six objects at z spe c < 2 is similar to that of the 20 remaining 
objects 



low specific SFRs. To resolve this issue it is necessary to 
combine our spectra with deep mid/far-infrared imaging 
to detect hidden star formation. 

It is interesting to compare our fraction of galaxies 
with strongly suppressed star formation to previous stud- 
ies 1 £^_Ojjt^j ; £^lt_isconsistent with the fraction found 
bv lLabbe et alJ ()2005[K as they identified three "red and 
dead" galaxies from a sa mple of 11 DRGs in the HDF- 
S, and with the study of IR eddv et al.l l)2006t) . who find 
that seven out of 24 DRGs at 1.5 < z < 2.6 are not de- 
tected at 24 /im, and have an average low specific SFR 
of 0.05 Gyr -1 . However, our fract ion is signif i cantly 
high er than the fractions found by IDaddi et alJ (2004) 
and IPap ovich e t al.l (J2006). Using the BzK criterion, 
IDaddi et alJ l)2004fr identifv all their 11 if -selected galax- 
ies at 2.0 < z < 2.7 as star-forming galaxies. The dif- 
ference may be partly explained by different definitions, 
as ~4 of our 9 galaxies with suppressed star formation 
galaxies would have been identifi ed as star-forming galax- 
ies by the BzK criterion. IPapovich et alJ l)2006|) finds 
that ~10% of 153 DRGs at 1.5 < z p hot < 3 show no signs 
of current star formation. Again, the selection criteria 
could play a role; no ne of our galaxies wo uld have been 
classified as "dead" bv lPapovich et all l)2006[) , as they ap- 
ply the following criteria: age > IGyr , E(B — V) < 0.1, 
age/T > 3 and no X-ray or 24/iin detection. The sam- 
ple selection could also be a factor, as the contribution 
of dusty star-forming galaxies to the DRG population is 
expected to be higher at z < 2. Further more, the CDF- 
S field - i n which both the studies by IPapovich et alJ 
( 2006) a nd IDaddi et al.l (12004ft wer e performed - may be 
atypical {van Dok kum et al.ll2006|) . We stress that our 
study is the first that is based on spectroscopic redshifts, 
and that this might also account for differences in the 
obtained fractions. This will be explored in a future pa- 
per. 

We note that most current CDM galaxy formation 
models fail to produce the high fr action of red galaxies a t 
the massive end as found bv Ivan Dokkum et alJ |2006). 
One easy way to solve this is by allowing more dust in 
the galaxies, but we have shown here that a large fraction 
of galaxies has low specific SFRs, and these are gener- 
ally a bsent in these models at z > 2 (e.g.. i Somervilld 
12001 iNagamine et al.ll2005l iKang et alJl2006l) . Our re- 
sults may indicate that the suppression of the star for- 
mation in the most massive galaxies occurs at higher 
redshift than has been predicted by current models. In 
this context, it is interesting to note that not all of the 
line-emission in our sample may be due to star forma- 
tion: from our low-resolution GNIRS spectra it appears 
that several of the emission-line galaxies with low spe- 
cific SFRs exhibit high [Nil] /Ha ratios, possibly indi- 
cating AGN. High resolution spectra of the line-emitting 
objects in our sample will be discussed in a future paper 
(M. Kriek et al., in preparation). 
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